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ABSTRACT 

We have obtained HST WFPC2 observations of three very luminous but mor- 
phologically different giant H II regions (GHRs) in M101, NGC5461, NGC5462 
and NGC5471, in order to study cluster formation in GHRs. Images obtained 
in the F547M and F675W bands are used to identify cluster candidates and for 
photometric measurements, and images in the F656N band are used to show ion- 
ized interstellar gas. The measured colors and magnitudes are compared with the 
evolutionary tracks generated by the Starburst99 and Bruzual & Chariot popula- 
tion synthesis models to determine the ages and masses of the cluster candidates 
that are more luminous than Mps47M = —9.0. The brightest clusters detected in 
the PC images are measured and found to have effective radii of 0.7-2.9 pc. 

NGC 5461 is dominated by a very luminous core, and has been suggested to 
host a super-star cluster (SSC). Our observations show that it contains three 
R136-class clusters superposed on a bright stellar background in a small region. 
This tight group of clusters may dynamically evolve into an SSC in the future, 
and may appear unresolved and be identified as an SSC at large distances, but 
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at present NGC5461 contains no SSCs. NGC5462 consists of loosely distributed 
H II regions and clusters without a prominent core. It has the largest number 
of cluster candidates among the three GHRs studied, but most of them are faint 
and older than 10 Myr. NGC 5471 has multiple bright H II regions, and contains 
a large number of faint clusters younger than 5 Myr. Two of the clusters in 
NGC 5471 are older than R136, but just as luminous; they may be the most 
massive clusters in the three GHRs studied. 

The fraction of stars formed in massive clusters has been estimated from the 
clusters' contribution to the total stellar continuum emission and from a com- 
parison between the ionizing power of the clusters and the ionizing requirement 
of the associated H II regions. Both estimates show that < 50% of massive stars 
are formed in massive clusters; consequently, the Ha luminosity of an H II region 
does not provide a sufficient condition for the existence of SSCs. The cluster 
luminosity functions (LFs) of the three GHRs show different slopes. NGC 5462 
has the steepest cluster LF and the most loosely distributed interstellar gas, 
qualitatively consistent with the hypothesis that massive clusters are formed in 
high-pressure interstellar environments. The combined cluster LF of the three 
GHRs has a slope similar to the universal cluster LFs seen in starburst galaxies 
and non-starburst spiral galaxies. 

Subject headings: galaxies: individual (M101) — galaxies: star clusters — H II 
regions — stars: formation 

1. Introduction 

Giant H II regions (GHRs) are sites of intense massive star formation. Their Ha lumi- 
nosities, 10 39 -10 41 ergs s _1 (Kennicutt 1984), require an ionizing power equivalent to that of 
24-2400 05V stars (Schaerer & de Koter 1997). With such high concentrations of massive 
stars, GHRs provide an excellent laboratory to study the modes of massive star formation, 
and in particular to probe whether they are birthplaces of globular clusters (Kennicutt & 
Chu 1988). 

In the two nearest GHRs, 30 Dor in the Large Magellanic Cloud (LMC) and NGC 604 
in M33, where stellar contents are well resolved, two distinct types of stellar groupings have 
been observed: 30 Dor is dominated by one central massive cluster R136 (Hunter et al. 1995; 
Walborn & Blades 1997), while NGC 604 contains multiple OB associations spreading over 
a large area (Hunter et al. 1996). Evidently, not all GHRs contain massive compact clusters; 
what physical environments give rise to the various cluster morphologies is currently under 
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investigation. 

One obvious way to elucidate this issue is to carry out detailed examinations of relatively 
nearby clusters and their environments. Maiz-Apellaniz (2001) studied 27 nearby (< 5 Mpc) 
clusters of varying morphological types, primarily classifying clusters based on their core and 
halo sizes. He suggests that compactness of clusters is predominantly related to the central 
density of the progenitor giant molecular cloud, i.e., extremely high pressure environments 
may be required to form massive compact clusters. However, this scenario has not been 
observationally tested; we do not know the pressures and densities of the giant molecular 
cloud in which optically visible clusters were formed. Examining clusters in a range of 
present-day environments may help us gain insight into their properties and relationship to 
their natal interstellar medium. 

Of all massive compact clusters, the most impressive ones are the super-star clusters 
(SSCs) with masses of 10 5 -10 6 M . SSCs are frequently observed in galaxy mergers and 
starburst galaxies (Whitmore 2003, and references therein), and they are believed to be 
preferentially formed in high-pressure interstellar conditions (e.g., Elmegreen & Efremov 
1997). However, some GHRs in non-interacting, late-type spiral galaxies may also host 
SSCs, especially those GHRs that are several times as luminous as 30 Dor and require ionizing 
powers rivaling those of young SSCs (e.g., Luridiana & Peimbert 2001). It is thus intriguing to 
examine the cluster content of such GHRs and investigate whether these relatively quiescent 
environments can also produce SSCs. 

The giant spiral galaxy M101 contains a large number of very luminous GHRs whose 
stellar content can be resolved and studied with Hubble Space Telescope (HST) images. We 
have therefore obtained HST observations of three M101 GHRs with different morphologies 
and galactic locations: NGC5461, NGC5462, and NGC5471. The locations of these GHRs 
in M101 are marked on a Second Palomar Observatory Sky Survey (POSS-II) red image in 
Figure 1. The properties of these GHRs are summarized in Table 1; for comparison, 30 Dor 
is also included in this table. We have used the HST continuum and Ha images of these 
three GHRs to carry out a detailed photometric study of their clusters. This paper reports 
our observations (§2) and methodology (§3), describes the cluster content in each GHR (§4), 
discusses cluster properties and their implications in studying massive star formation and 
cluster formation (§5), and summarizes our results (§6). 



-4- 



2. Observations and Data Reduction 

The HST WFPC2 images of the GHRs NGC5461, NGC5462, and NGC5471 were 
obtained for the Cycle 6 program GO-6829. The observations were made through the con- 
tinuum filters F547M (Strdmgren y) and F675W (WFPC2 R), and the Ha filter F656N (for 
filter characteristics, see Biretta et al. 1996). Multiple exposures in each filter were made 
with a GHR centered on a Wide Field Camera (WFC) for photometric measurements. Two 
short exposures in F547M with the GHR centered on the Planetary Camera (PC) were also 
made for cluster size measurements. The observations are listed in Table 2. 

We received the HST pipeline processed WFPC2 images and then reduced them fur- 
ther with the IRAF and STSDAS routines. All images were corrected for the intensity- 
and position-dependent charge transfer efficiency (CTE) by applying a linear ramp with a 
correction factor chosen according to the average counts of the sky background (Holtzman 
et al. 1995). Images in each filter were then combined to remove cosmic rays and to produce 
a total-exposure map. To better illustrate the spatial correlation between the stars/clusters 
and the ionized gas, we have produced color images of the three GHRs using a customized 
IDL routine. The individual F547M, F675W, and Ha images were mapped to the colors 
blue, green, and red, respectively. These images were transformed to a logarithmic scale, 
and the maximum and minimum flux values for each filter were chosen in order to maximize 
the dynamic range of the image while also creating a relatively black background. The color 
images of NGC5461, NGC5462, and NGC5471 are shown in Figures 2-4, where the ionized 
gas appears red and most stars blue. The individual F547M, F675W, and Ha images of 
NGC5461, NGC5462 and NGC5471 are presented in Figures 5-7, respectively. 

Aperture photometry was carried out using the IRAF task APPHOT for the F547M and 
F675W images. Owing to the small number of identifiable candidate sources and the complex 
blending and irregular background in some regions, we manually selected compact sources 
with obvious peaks as candidate clusters in the three GHRs. The candidate clusters in the 
three GHRs are marked in Figures 5-7. The apparent magnitudes, mp547M and mp675W> 
were measured with the WFC images using a source aperture of radius 2 pixels (0'/2) and 
an annular background aperture of radii 6-11 pixels. For clusters with neighboring clusters 
within < 0"3, such as #8, #9 and #10 in NGC5461 and #1 and #2 in NGC5471, the 
photometry was measured with a O'.'15-radius source aperture using both the WFC and 
PC images. The corrections from the (y/15-radius aperture to the 0V2-radius aperture are 
determined by measuring isolated resolved and unresolved sources in each image, and are in 
the range of —0.2 to —0.3 mag. The fluxes of these clusters will be over-estimated from the 
WFC images due to the inclusion of neighbor's light, and the errors are larger for the fainter 
clusters. For example, the error in the mp675w of NGC 5461-9 may be as large as —0.4 mag, 
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which is estimated by comparing its 7BF547M measured from the PC and WFC images. 

We have derived the magnitudes in the VEGAMAG system. The errors in m F547M 
given by APPHOT are ~ 0.01 mag for m F547M < 20, and rise to 0.02-0.03 mag for m F5 4 7M = 
20-21. The errors in m F675W are generally larger, with most of them < 0.02 mag but 
some as high as 0.04 mag. These formal errors are derived from the flux variations in the 
background annulus used in the photometric measurements. In regions with bright irregular 
backgrounds, using different background apertures may produce different photometric results 
and the uncertainties in photometry will be larger than the formal errors given by APPHOT. 
We have taken these uncertainties into account and estimated realistic errors. 

The bright irregular sky background is attributed to both stars and nebulosity. The 
extended distribution of unresolved stars, similar to the "star clouds" defined by Lucke & 
Hodge (1970) for OB associations in the LMC, raises the background in both the F547M and 
F675W images. We find that using different annular sky backgrounds results in uncertainties 
of ~ 0.03-0.05 mag in m F54 7M for sources near modest star clouds, ~ 0.1 mag for sources 
surrounded by bright star clouds, and up to ~ 0.2 mag for faint sources near bright star 
clouds. The stellar background does not affect the (m F5 47M — ^F675w) color as much because 
the variations in m F54 7M and m F675 w are correlated. The bright nebular background, on 
the other hand, contributes uncertainties only to m F 675w and hence affects the color. The 
uncertainties in the (m F 547M — ^F675w) color are ~ 0.05-0.06 mag for most sources near 
nebulosities, and up to ~ 0.16 mag for faint stellar sources near bright nebulosity. 

To reduce the uncertainties introduced by nebular contamination in the F675W im- 
ages, we have produced Ha-free F675W images by subtracting scaled Ha images from the 
F675W images. The Ha-subtracted F675W images of the three GHRs are also presented 
in Figures 5-7. Aperture photometry has been carried out for the Ha-subtracted F675W 
images, and the apparent magnitude is designated as mp675W'- The uncertainties in the 
(m F54 7M — ^F675W') color with different annular sky backgrounds are reduced to ~ 0.02-0.05 
mag for most candidate cluster sources near nebulosities. 

Our O'.'2-radius source aperture does not include all the light from a cluster. The correc- 
tion from a C.^-radius aperture to a O'.'5-radius aperture, which includes ~ 95% of the light 
of a point source (Holtzman et al. 1995), is ~ —0.20 ± 0.06 mag for unresolved sources and 
larger for resolved sources. As the cluster candidates are resolved to different extents, the 
aperture corrections are in the range of —0.2 to —0.3 mag but difficult to determine exactly. 
We have chosen not to apply aperture corrections; therefore, our photometric measurements 
are systematically fainter by 0.2-0.3 mag, but the analysis and conclusions of this paper are 
not sensitive to such small offsets that incur on the magnitudes and colors of the clusters. 
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The photometric results of NGC5461, NGC5462 and NGC5471 are presented in Ta- 
bles 3-5, and plotted in the color-magnitude diagram (CMD) of MF547M versus (MF547M — 
Mp675w) i n Figures 8-10, respectively. For clusters with close neighbors, their Mps47M mea- 
sured from the PC images and their (Mps47M — ^F675W') measured from the WFC images 
are used in the CMD. These absolute magnitudes are derived using a distance modulus of 
(m — M) = 29.3 (Stetson et al. 1998). The Galactic foreground extinction toward M101, 
E(B — V) = 0.01 (Schlegel, Finkbeiner, & Davis 1998), is corrected, although its effect is 
negligible. The internal extinction from M101 is not individually corrected for, given that it 
is highly variable and the measurements are only available for certain parts of the GHRs. 



3. Methodology 

The observed magnitudes and colors of the clusters can be used to determine their ages 
and masses through comparisons with those predicted by population synthesis models (e.g., 
Elson & Fall 1985; Bruzual & Chariot 1993, 2003; Leitherer et al. 1999). Below we describe 
the synthetic photometry derived from models and how we use it to estimate the properties 
of clusters. We have also used the Larsen (1999) method to measure the sizes of clusters from 
their surface brightness profiles. The procedures of cluster size measurements are outlined 
at the end of the section. 



3.1. Synthetic Photometry 

We have used the Starburst99 models (Leitherer et al. 1999) and the Bruzual & Chariot 
(2003, hereafter BC03) models to generate synthetic photometry for comparison with obser- 
vations of our cluster candidates in GHRs. We have adopted a Salpeter initial mass function 
(IMF) with lower and upper mass limits of 1 M Q and 100 M , which are commonly used in 
population synthesis models for star-forming and starburst regions. The luminosity, colors, 
and evolution of a cluster depend on its metallicity. To select appropriate models, we have 
used the observed oxygen abundances of the GHRs to assess their metallicities, because 
clusters and their surrounding GHRs are expected to have the same abundances and the 
oxygen abundances are well determined. Oxygen abundances of NGC5461 and NGC5462, 
relative to the solar value, have been measured to be 0.6-0.9, while that of NGC5471 is ~ 
0.25 (e.g., Evans 1986; Scowen, Dufour, & Hester 1992; Pilyugin 2001; Luridiana et al. 2002). 
Therefore, we adopt the 1 Z Q model for NGC5461 and NGC5462, and the 0.2 Z Q model for 
NGC5471. 
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The F547M and F675W niters we used are not included in the default filter systems of 
Starburst99 or BC03 for which synthetic photometry is readily available; thus, customized 
procedures are needed to derive synthetic Mp547M and Mf675w- As the first step, we use 
Starburst99 Version 4.0 to generate integrated stellar spectra for a simple stellar population 
(SSP; i.e., a single-age and single-abundance group of stars) from ages of to 30 Myr at 1 
Myr intervals and from 30 to 150 Myr at 3 Myr intervals. This step is not necessary for the 
BC03 models, as integrated spectra for an SSP are available for most of these age intervals. 
These model spectra are those without nebular line and continuum emission because the 
clusters are generally well resolved from the superposed extended nebular emission and the 
background-subtraction in APPHOT adequately removes the extended nebular emission. 
The synthetic spectra from Starburst99 and BC03 are then convolved with filter transmission 
curves, using the IRAF/STSDAS task CALCPHOT, to calculate the synthetic Mp 547 M and 
Mf675W- We have produced synthetic photometry for SSPs with metallicities of 0.2 Z & and 
1 Zq, and generated evolutionary tracks in the CMDs in Figures 8-10 for comparisons with 
observations of NGC5461, NGC5462, and NGC5471, respectively. The differences in the 
two sets of evolutionary tracks reflect the differences between the Geneva and Padova stellar 
evolution models used by Starburst99 and BC03, respectively. However, the effects of these 
differences are small compared to the uncertainties in the cluster mass estimates. 

We use the R136 cluster at the core of 30 Dor as a reference point, because it is an 
archetypical populous blue cluster and possibly a young globular cluster. The R136 cluster is 
~ 3 Myr old (Hunter et al. 1995; Walborn & Blades 1997). Its spatially resolved photometry 
in the V band has been measured and the absolute visual magnitude within a radius of 7 pc is 
My = —11.1 (Moffat, Seggewiss, & Shara 1985). This V band magnitude is adopted directly 
because its central wavelength is similar to that of F547M and the 7-pc-radius aperture 
matches that used in the photometric measurements of M101 clusters. The (Mps47M — 
Mf675w) color of R136 is not available, so we use the synthetic color generated by Starburst99 
for a 3 Myr old cluster with Z = 0.2-0.4 Z@. As no extinction correction has been applied 
to the M101 clusters, we have reddened the synthetic color of R136 with its visual extinction 
A v = 1.2 (Moffat et al. 1985) and marked both the dereddened and reddened R136 in the 
CMDs in Figures 8-10. 



3.2. Assessing Masses and Ages of Clusters 

The mass and age of a cluster can be assessed by comparing its magnitudes and colors to 
model predictions if photometric data are available in three passbands. For a young cluster, 
it is important to include a U band or a B band because the spectral energy distribution of 
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young massive stars peaks in the ultraviolet wavelengths. Unfortunately, only F547M and 
F675W photometry is available for the clusters in M101 GHRs, and these two bands are 
not as sensitive to the young massive stars as the U and B bands. We cannot determine 
unambiguously the cluster masses and ages by comparing the photometric measurements 
with the evolutionary tracks of SSPs in the CMDs. However, we may use the interstellar 
environment as an independent diagnostic of the cluster age, as the interstellar medium 
around a cluster evolves as a result of stellar energy feedback. 

At ages < 5 Myr, a cluster has the highest ionizing power, and hence will be in a dense, 
luminous H II region. At ages 5-10 Myr, the fast stellar winds and supernova explosions from 
a cluster have swept up the ambient ISM into a supershell with a visible cavity around the 
cluster. At ages > 10 Myr, a cluster loses its ionizing power and has dispersed its ambient gas, 
so it will be surrounded only by diffuse gas with low surface brightness. We have compared 
the Ha images with the continuum images to examine the interstellar environment of the 
clusters and to assess the approximate ages of the clusters. In Figures 8-10, we mark circles 
around the clusters that are coincident with compact, luminous H II regions, and dashed 
circles around the clusters that are in supershells, indicating that their ages are < 5 Myr 
and 5-10 Myr, respectively. The unmarked clusters, not surrounded by bright Ha emission, 
are older than 10 Myr, but their exact ages are poorly constrained. 

With a rough estimate of the cluster age, it is then possible to compare the location 
of a cluster in the CMD with the synthetic evolutionary tracks of clusters to determine the 
cluster mass. The photometric measurements of the clusters have not been corrected for the 
extinction within M101, thus using these reddened magnitudes would underestimate cluster 
masses. To illustrate the effect of extinction and to make a rough correction, we take the 
visual extinctions of the GHRs determined from their Balmer decrements by Kennicutt & 
Garnett (1996), and plot the corresponding reddening vectors in Figures 8-10. We have 
adopted these nebular extinctions and made reddening-corrected estimates of masses for 
clusters more luminous than Mp547M — — 9. The age and mass estimates of these luminous 
clusters are given in Table 6. We do not attempt to estimate masses for fainter clusters 
because these luminosities overlap those of single supergiants (Humphreys & Davidson 1979). 
Furthermore, many faint clusters are not surrounded by bright nebulosity, indicating poorly 
constrained ages at >10 Myr, so their mass estimates would be highly uncertain. 

We have used the R136 cluster to estimate the uncertainties in our cluster mass esti- 
mates. From the extinction-corrected location of R136 in the CMD, we estimate a mass of 
~ 2 x 10 4 M . The mass of R136 has been derived from its resolved stellar content to be 
2.2 xlO 4 M (Hunter et al. 1995) by summing the masses of stars > 2.8 M (mass cutoff 
limited by completeness) within a 4.7-pc radius. Note however that our estimate of mass is 
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based on the luminosity of R136 within a 7-pc radius and a minimum stellar mass of 1 M . 
Our mass estimate of R136 using the 4. 7-pc radius aperture (Moffat et al. 1985) and the 2.8 
M lower mass limit is 1.4 xlO 4 M , about 40% lower than that derived from the resolved 
stellar content. Therefore, the uncertainties in our cluster mass estimates are at least 40%. 

3.3. Assessing Cluster Sizes 

Some of our clusters appear resolved in the PC images, so it is possible to determine 
their sizes. The size of a cluster can be described by its effective radius, R c g, the radius 
that encircles half of the cluster light. The R c g of a cluster can be estimated with the 
routine ISHAPE developed by Larsen (1999). In this routine, the surface brightness profile 
of a cluster is modeled by an analytic function and convolved with a point spread function 
(PSF) calculated with the TINY TIM Version 6.0 (Krist 1995) for the cluster's position on 
the PC chip. The PSF-convolved model profile is then compared with the observed cluster 
profile. The best-fit model, judged by the x 2 statistics, gives the full width at half maximum 
(FWHM) of the analytic function, which is then used to determine R e s. 

We have used the two most common analytic functions, the Gaussian and the King 
(1962) profiles, to model the clusters. The King profile contains a concentration parameter 
c = log(r t /r c ), where r t is the tidal radius and r c is the core radius. Typically c is within 
the range of 1.0-2.0 for globular clusters in the Galaxy (Harris 1996) and young rich clusters 
in the LMC (Elson, Fall, & Freeman 1987). We have experimented with different values 
of c within this range in the model fits, and found that for a cluster detected with S/N > 
15, the best-fit R e s is insensitive to the concentration parameter c or even the form of the 
analytic function. For a bright cluster with adequate S/N, the R c g is estimated using both 
the Gaussian and the King profiles, and the average of the two estimates is adopted and 
given in Table 6. 

All but one of the cluster sizes we measured are in the range of i? e fr = 0'. / 02-0 / .'09, 
corresponding to 0.7-2.9 pc. (Note that R c s can be smaller than the pixel size of 0'/0455, 
because the PSF effects have been considered and removed in the profile fitting.) The only 
exception is NGC 5471-9, whose R e{i is only 0"005; as we discuss later in §4.3, NGC 5471-9 
is most likely a luminous A-F supergiant, instead of a cluster. The sizes of these M101 
clusters are within the range of clusters in the Galaxy and nearby galaxies. For example, the 
globular clusters in the Galaxy have _R c fr ~ 1-5 pc, with a median of ~3 pc (Harris 1996); 
the compact young cluster R136 in the LMC has i? e fr ~ 0.9 pc (Mackey & Gilmore 2003) 1 ; 



The surface brightness profile of R136 indicates a compact, dominant component on top of a broad, 
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compact young massive clusters in nearby starburst galaxies where stars are not resolved 
have R eS ~ 2-4 pc (Meurer et al. 1995). Comparisons between the M101 clusters and R136 
will be discussed in more detail in §4. 

4. Clusters in Three Luminous GHRs in M101 

Below we describe the spatial distribution, ages, masses, and sizes of the clusters in 
NGC5461, NGC5462, and NGC5471. The extinctions of individual H II regions in these 
GHRs are taken from Kennicutt & Garnett (1996). 

4.1. Clusters in NGC 5461 

The GHR NGC 5461 has been loosely defined to be the H II complex extending over 
a 66" x 26" region in ground-based Ha images (Israel, Goss, & Allen 1975). Considering 
that this area corresponds to a linear size of ~ 2.3 kpc x 0.9 kpc, it is unlikely that the 
entire region is associated with one coherent star formation event. Indeed, 12 H II regions 
have been identified within NGC 5461 by Hodge et al. (1990). Our WFPC2 Ha image shows 
that NGC 5461 contains two regions that would have been individually identified as GHRs 
if they were in the Local Group: H1105 and H 1098 (marked in Figure 5e; designation from 
Hodge et al. 1990). H1105 is 3 times as luminous as 30 Dor, and H1098 is as luminous as 
NGC 604, or 1/3 as luminous as 30 Dor. The 10 fainter H II regions are distributed roughly 
along the axis connecting H 1105 and H 1098 with a higher concentration toward H 1105. We 
define the "main body" of NGC 5461 to be the region containing H1105 and H 1098 and 
their vicinity, as in the field-of-view of Figure 5. 

The relationship between the stars/clusters and the H II regions is clearly illustrated in 
the color image of NGC 5461 (Figure 2). In addition to the bright H II regions, NGC 5461 
also has nebular filaments, loops, and well-defined shells with stars/clusters underneath 
these interstellar structures. A total of 12 cluster candidates are identified in the main body 
of NGC 5461; they are listed in Table 3 and marked in Figure 5f. Only six clusters have 
M F547 m < -9: five in H1105 (#6, #8, #9, #10, and #11) and one in H 1098 (#1). A 
careful inspection of their immediate surroundings shows that all six clusters are superposed 
on bright H II regions, indicating ages of < 5 Myr. To estimate the masses of these young 



shallow component. The R e g is estimated using its core radius of 0.32 pc and a King profile with c = 1.5 to 
approximate the compact component. 
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luminous clusters, we adopt the visual extinctions of A v — 1.7 ± 0.4 and 0.8 ±0.1 of H 1098 
and H 1105, apply the respective extinction correction, and compare the dereddened cluster 
positions in the CMD in Figure 8 with the evolutionary tracks generated by Starburst99 
and BC03 for different SSP masses. Four of the clusters, #1, #6, #9, and #10, show 
dereddened colors consistent with SSPs at ages < 5 Myr; thus their masses can be estimated 
in a straightforward manner. The two remaining clusters, on the other hand, have dereddened 
colors consistent with SSPs at ages greater than 6 Myr. The red colors of clusters #8 and #11 
can be caused by large local extinction excesses or stochastic color deviation for low-mass 
SSPs. We consider the latter more likely, i.e., the cluster contains or is projected near a red 
supergiant and the cluster color is thus confused. For example, cluster #11 may consist of 
(or be projected towards) a K0 Ia-0 supergiant with My = —9.4 (Humphreys 1978) and (V- 
R)o = 0.76 (Johnson 1966) and a young cluster with Mps47M — —9.2. The masses of clusters 
#8 and #11 are determined with the ad hoc assumption of a contaminating red supergiant, 
which reduces the luminosity of the cluster and lowers the mass estimate accordingly. The 
mass estimates of the six brightest clusters in NGC5461 are mostly in the range of 1-3 xlO 4 
M Q (see Table 6). These masses are comparable to that of R136, ~2x 10 4 M Q . 

The PC images of NGC5461 are used to determine the cluster sizes, but only three 
clusters, #6, #8, and #10, are detected with S/N >15 for reliable size measurements. The 
R c s estimated for these three clusters are 0.8, 0.7, and 2.1 pc, respectively. While the sizes 
of clusters #6 and #8 are comparable to that of R136, R c g = 0.9 pc, cluster #10 is more 
extended and shows visible departure from spherical symmetry in the PC image in Figure 11. 
The morphology of cluster #10 suggests that it may be a composite of two clusters with the 
southwest object brighter than the northeast object. 

Among the three GHRs we studied in M101, NGC5461 is particularly interesting be- 
cause it is one of the most luminous GHRs in galaxies within 10 Mpc (Kennicutt 1984). 
Furthermore, the core of NGC5461 (i.e., H1105) has a remarkably high surface brightness 
with a peak emission measure of 4.4 x 10 5 cm" 6 pc, comparable to those of the most ac- 
tive starburst regions. The Ha luminosity of H1105 implies an ionizing flux rivaling those 
of SSCs (Kennicutt & Chu 1988; Luridiana & Peimbert 2001); thus, NGC5461 has been 
considered the most promising site in M101 where SSCs might be found. However, our 
analysis shows that H 1105 contains five R136-class clusters, which are by no means in the 
same league as the SSCs with masses ~ 10 5 -10 6 M commonly found in starburst galaxies 
or mergers (e.g., O'Connell, Gallagher, & Hunter 1994; Whitmore et al. 1999). The core of 
NGC 5461 is nevertheless striking in its high concentration of stars in a small volume - three 
clusters with a total mass of 6 x 10 4 M in a region of ~ 32 pc (0'/9) across. It is possible 
that these clusters are subclusters that will dynamically interact and merge into a cluster 
that has a mass more typical for SSCs. As shown in the numerical simulations of Bonnell, 
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Bate, & Vine (2003), the hierarchical fragmentation of giant molecular clouds naturally leads 
to the formation of subclusters that can merge to form the final stellar cluster. If H1105 
were projected to a distance similar to that of the Antennae galaxies (~ 20 Mpc, Whitmore 
et al. 1999), it would imitate a single cluster as shown in Figure lib, and the combined 
light of the clusters and the bright stellar background in a O^-radius aperture would have 
A^F547M = —13.0, corresponding to a mass of ~ 10 5 M . These results suggest that some 
young SSCs previously identified at distances of > 20 Mpc may be tight groups of R136-class 
clusters as seen in the core of NGC5461. 



4.2. Clusters in NGC 5462 

The GHR NGC 5462 corresponds to a large H II complex with a dimension of 90" x 34", 
or 3.2 kpc x 1.2 kpc, in ground-based Ha images (Israel et al. 1975). Thirty-three H II regions 
have been identified within NGC 5462 (Hodge et al. 1990), but none are comparable to 30 
Dor. The overall morphology of NGC 5462 consists of a few bright H II regions distributed 
along an axis from northeast to southwest and fainter filaments and loops extending outwards 
from this axis. The two brightest H II regions, H 1170 and H 1176, are each only comparable 
to NGC 604, or 1/3 as luminous as 30 Dor; the others are much fainter. The distribution of 
star formation in NGC 5462 is apparently not as concentrated as in NGC 5461. We define 
the "main body" of NGC 5462 to be the region containing H1176, H1170, H1159 and their 
vicinity, as in the field-of-view of Figure 6. 

The color image of NGC 5462 (Figure 3) shows a distinct offset between the ionized 
gas and concentrations of stars, suggesting that the star formation has proceeded from the 
southeast to northwest. A total of 25 loosely distributed cluster candidates are identified in 
NGC 5462 (see Table 4 and Figure 6f). Most of these clusters are faint; only three clusters, 
#6, #18, and #23, have M F 547M < —9 and are analyzed for their masses. While cluster 
#18 is superposed on a bright H II region H 1176, indicating an age of < 5 Myr, clusters 
#6 and #23 are not associated with any H II regions or supershells, indicating ages > 10 
Myr. For cluster #18 the visual extinction Ay = 0.9 ± 0.4 of the surrounding H II region 
H 1176 is adopted, and for the other two clusters the visual extinctions of their nearest H II 
regions are adopted, i.e., Ay = 0.6 ± 0.2 of H 1159 for cluster #6 and Ay = 0.9 ± 0.4 of 
H1176 for cluster #23. We apply the respective extinction correction to each of the three 
bright clusters and compare their dereddened cluster positions in the CMD in Figure 9 with 
the evolutionary tracks. The young cluster #18 shows a dereddened color consistent with 
SSPs at ages of < 5 Myr and the estimated mass is < 1 x 10 4 M . The two older clusters 
#6 and #23 show dereddened colors consistent with SSPs at ages of > 10 Myr; however, as 
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their ages are poorly constrained, we only obtain their lower mass limits by assuming cluster 
ages of ~ 10 Myr. We note that unless the two clusters have ages > 30 Myr, their masses 
would be within a factor of 1.5 of the lower mass limits. The mass estimates of the three 
brightest clusters in NGC5462 are in the range of 1-2 xlO 4 M (see Table 6), comparable 
to the mass of R136. 

In the PC images of NGC5462, only cluster #6 is detected with S/N > 15 for size 
measurements. The i? e fr estimated for this cluster is 2.3 pc, more extended than that of 
R136. As shown in Figure 12, cluster ^6 has an asymmetric morphology elongated along 
the northwest and southeast direction, indicating a complex structure. 

NGC5462 has the lowest Ha surface brightness among the three GHRs we studied 
in M101. It has a larger number of clusters than the other two GHRs, but only three are 
R136-class clusters. Furthermore, the clusters do not show obvious spatial concentrations, in 
sharp contrast to those seen in NGC 5461. The combination of the low Ha surface brightness 
and sparse distribution of small-mass clusters suggests that the star formation and cluster 
formation is more spread-out and modest in NGC 5462. 

4.3. Clusters in NGC 5471 

The GHR NGC 5471 extends over a diameter of ~ 17", or ~ 600 pc. Ground-based 
images of NGC 5471 show five bright knots, which are designated as A, B, C, D, and E 
components by Skillman (1985) and have been called NGC5471A-E, respectively. Our color 
image of NGC 5471 (Figure 4) shows that the A-, B-, C-, and E-components display bright 
H II regions centered on clusters. However, the D-component displays an offset between the 
H II region and the clusters, which are located to the north and the east sides of a dark 
cloud, respectively; it is uncertain whether the clusters and the H II region are physically 
associated. The A-component is as luminous as 30 Dor, and the B-, C-, and E-components 
are comparable to or fainter than NGC 604. 

A total of 19 cluster candidates are identified in NGC 5471; they are listed in Table 5 
and marked in Figure 7f. Most of the clusters reside within the A-E components, with 
the highest concentration located in the A-component and its western extension. The eight 
clusters with Mp547M < — 9 have been analyzed. We have adopted the visual extinction of 
each component, applied the respective extinction correction to the eight brightest clusters, 
and compared the dereddened cluster positions in the CMD with the evolutionary tracks (see 
Figure 10). The comparisons are problematic because most of the dereddened cluster colors 
do not agree with those expected from the evolutionary tracks. While the disagreements are 
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partially attributed to errors in the extinction and photometric measurements, the dominant 
cause of the disagreements is probably uncertainties in the stellar evolution models at low 
metallicities. It is known that at Z < 0.2 Z Q stellar evolution models cannot reproduce 
the observed luminosities and colors of red supergiants or the number ratios of blue to red 
supergiants (Mayya 1997; Origlia et al. 1999; Leitherer et al. 1999). These uncertainties 
directly affect the luminosities and colors of SSPs, particularly at ages around 7 to 14 Myr 
when red supergiants are significant contributors of the total light. Given these uncertainties, 
we can make only order-of- magnitude mass estimates for the clusters in NGC5471. 

Four of the clusters we analyzed (#2, #7, #12, and #16) are superposed on bright H II 
regions, suggesting that they are < 5 Myr old. Clusters #7, #12, and #16 are only ~ 0.1 mag 
bluer or redder than those expected for young clusters; therefore, we disregard these color 
differences and use only the M F5 4 7M to estimate their cluster masses. Cluster #2, on the other 
hand, is ~ 0.6 mag redder than the color expected for its young age, and this discrepancy 
is larger than the known errors in photometry or stellar evolution models. We suggest that 
this red color excess is likely attributed to a contaminating post-outburst luminous blue 
variable (LBV) because cluster #2 is located in the C-component where high-velocity (> 
1000 km s" 1 ), [N II]-bright nebular emission similar to that of i] Car's ejecta nebula has 
been reported (Castaneda, Vilchez, & Copetti 1990). We assume that cluster #2 contains 
an LBV similar to rj Car, which has V = 6.22 and R = 4.90 at quiescent states (Mendoza 
1967) and can brighten up by 1-2 mag during outbursts or 3-5 mag during super-outbursts 
(Humphreys & Davidson 1994). These quiescent V and R magnitudes can be converted to 
^F547M = —8.5 and M F675W = —9.4 using the distance modulus of (m — M) = 12.79 and 
the extinction of Ay = 1.92 for rj Car's host cluster Trumpler 16 (DeGioia-Eastwood et al. 
2001). The remaining members of cluster #2 would have (Mp547M — ^F675w) — 0.44, which is 
still too red for a < 5 Myr old cluster. However, if the hypothesized LBV is 0.4 mag brighter 
(during or after an outburst), the rest of cluster #2 would have (M F 547m — ^F675w) — —0.1 
as expected for a young cluster and M F 547M = 8.2. These resultant color and magnitude are 
used to estimate the mass of cluster #2. 

The other four clusters we analyzed (#3, #4, #5, and #9) are not associated with 
H II regions or supershells, suggesting that their ages are > 10 Myr. However, as discussed 
in the next paragraph, #9 has a small intrinsic size which makes it more likely a luminous 
supergiant rather than a cluster. Among the remaining clusters, #4 shows a dereddened color 
consistent with SSPs at ages > 10 Myr and thus its lower mass limit is easily estimated. 
Clusters #3 and #5, on the other hand, have dereddened colors ~ 0.2 mag bluer than 
that of SSPs at ages > 10 Myr. Since their extinction corrections are already small, the 
disagreements are unlikely to be caused by the uncertainty in extinction measurements. It 
is most likely that the disagreements arise from the uncertainties in the modeled colors and 
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photometry, so we disregard the blue color excesses when estimating the cluster masses. 
The mass estimates for the eight bright clusters are in the range of ~ 0.5-2 xlO 4 M (see 
Table 6), approaching or comparable to R136. 

The cluster sizes are determined using the PC images of NGC5471. Four clusters, #4, 
#5, #9, and #16, are detected with S/N >15. The i? e fr estimated for these four clusters are 
2.9, 1.4, 0.2, and 1.1 pc, respectively. The small size of cluster #9 suggests that it is either 
an unresolved star or a post-core-collapse globular cluster (Harris 1996). Post-core-collapse 
globular clusters are at least 10 9 yr old (e.g., Heggie 1985) and thus unlikely to exist in GHRs. 
We consider it more likely that cluster #9 is a star with magnitudes and colors compatible 
to those of a luminous A-F supergiant (Moffat & Fitzgerald 1977; Humphreys, Massey, & 
Freedman 1990). For the three resolved clusters, #16 has a size comparable to that of R136, 
while #4 and #5 are more extended and show asymmetric, elongated morphologies in the 
PC image in Figure 12. 

NGC5471 has a large number of young clusters with ages < 5 Myr. However, the 
majority of these young clusters are faint with M F 547M > _ 9, which may be small clusters 
with masses of a few xlO 3 M or just luminous supergiants. Among the cluster candidates 
in NGC5471, #4 and #5 are the most massive ones since they are older than 10 Myr and 
still as luminous as the young R136 cluster, suggesting that their masses are higher that that 
of R136. 



5. Discussion 

5.1. Nature of Faint Cluster Candidates in GHRs 

A large number of cluster candidates have been identified in the three M101 GHRs, but 
the nature of the faintest objects is uncertain. It is possible that some of these faint cluster 
candidates consist of multiple OB associations as observed in the nearby GHR NGC604 
(Hunter et al. 1996) and some are simply luminous supergiants frequently seen near high 
concentrations of massive stars, such as 30 Dor (Walborn & Blades 1997). We have therefore 
simulated WFPC2 images of 30 Dor and NGC604 at a distance of 7.2 Mpc and searched 
for "clusters" using the same criteria as we did for cluster candidates in the M101 GHRs. 
The spurious clusters in 30 Dor and NGC 604 can be identified because their resolved stellar 
contents are known. The real and spurious clusters in these two GHRs can then be compared 
with the cluster candidates in M101 to better assess the nature of the latter. 

To simulate a WFPC2 image of 30 Dor at 7.2 Mpc, we have used a green continuum (A c 
= 5130 A, AA = 155 A) image from the Magellanic Clouds Emission-Line Survey (MCELS, 
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Smith & The MCELS Team 1999), and binned the data to 3.5 pc per pixel. The resultant 
image is displayed in Figure 13a. Using the same identification criteria for clusters in M101 
GHRs, the two brightest objects in 30 Dor will be selected as clusters: R136 and R131. 
While R136 is a bona fide cluster, R131 (= HD 269902) is an AO supergiant with V = 10.0 
(corresponding to My + Ay = —8.5) and hence a spurious cluster. The mis-identification of 
R131 as a cluster bolsters our choice of a luminosity cutoff of MF547M = —9.0 for cluster mass 
estimates (§3.2), as the fainter cluster candidates may be single supergiants. It is interesting 
to note that two other concentrations of stars in 30 Dor are not identified as clusters: the 
Hodge 301 cluster (Hodge 1988) and the OB association LH99 (Lucke & Hodge 1970). The 
Hodge 301 cluster, with an age of ~ 20-25 Myr and a mass of a few 10 3 M (Grebel & Chu 
2000), is too faint to meet our cluster identification criteria. LH99, on the other hand, is 
too distributed to mimic a cluster. 

A WFPC2 image of NGC604 at 7.2 Mpc is simulated with its archival WFPC2 F547M 
image. Adopting a distance of 0.84 Mpc to M33 (Freedman, Wilson, & Madore 1991), the 
data are binned to 3.5 pc per pixel, and the resultant image is displayed in Figure 13b. 
The OB associations in NGC604 appear as small concentrations on top of an irregular 
stellar background. The four brightest concentrations have Mps47M — —8.0 to —9.0, lumi- 
nous enough to meet our identification criteria for M101 clusters. However, most of the 
concentrations have irregular shapes, and the brightest one does not even have an obvious 
boundary. Therefore, the NGC604-type OB associations may mimic faint clusters at best, 
with luminosities rivaled by those of supergiant stars. 

For a direct comparison between M101 GHRs and the simulated 30 Dor and NGC604 
at 7.2 Mpc, Figure 13 displays their images in the same spatial and intensity scales 2 over a 
350 pc x 350 pc field-of-view. It is immediately clear that the brightest cluster candidates 
in M101 GHRs are more luminous than R136 and are most likely bona fide clusters. The 
nature of the cluster candidates fainter than M F 547M = —9.0, marked in Figure 13, is less 
obvious. Some faint cluster candidates may be blue supergiants because they have sharp 
images and appear isolated, and the blue color excludes the possibility of post-core-collapse 
clusters; examples of these include NGC 5461-7, NGC 5462-20, and perhaps NGC 5471-13 and 
14. Some faint cluster candidates may be OB associations because they appear extended 
without a sharp boundary; examples of these include NGC 5461-5, NGC 5462-13 and 15, and 
NGC 5471-10, 11, and 17. 



2 The green image of 30 Dor was taken in a similar but not identical wavelength band. The intensity scale 
of the 30 Dor image is selected to match that of the F547M images as much as possible so that objects with 
similar magnitudes appear similar in both F547M and green band images. 
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5.2. The Fraction of Massive Stars Formed in Clusters 

It has been suggested that massive stars form preferentially in associations and clusters 
(Stahler, Palla, & Ho 2000, and references therein). Our HST WFPC2 images of the M101 
GHRs show that several young R136-class clusters are superposed on discrete regions of 
unresolved stellar emission, e.g., clusters #8, #9, and #10 at the core of NGC5461 and 
cluster #16 at the core of NGC5471A. The unresolved stellar backgrounds are most likely 
star clouds that contain field stars and loosely assembled associations. The similarities in 
locations and colors suggest that the clusters and the background star clouds are formed 
from the same episode of star formation. It is then interesting to determine the fraction 
of massive stars that are formed in R136-class clusters in these regions. We have used two 
different methods to determine this fraction: one is based on the contribution of cluster 
light to the total light, and the other is based on a comparison between the ionization flux 
expected from the clusters and the ionizing flux required by the surrounding H II region. 

We have selected four regions for this analysis: two in NGC5461, one in NGC5462, 
and one in NGC5471. These regions are listed in Table 7, and their close-up F547M and 
Ha images are presented in Figure 14. The F547M images show that the clusters in these 
four regions are all superposed on discrete bright diffuse stellar backgrounds, and the Ha 
images show that all are at the cores of bright compact H II regions. We have measured the 
total light from these four regions in both F547M and F675W bands using the apertures 
marked on the F547M images in Figure 14 and described in column 3 of Table 7. The 
background, determined from the median of an annular region outside the H II region, has 
been subtracted, although it contributes to only 1-2% of the total light. For the clusters, we 
have applied aperture corrections of ~ —0.2 mag to our photometric measurements made 
with a O'.'2-radius APPHOT aperture to account for the missing light. The cluster-to-total 
light ratios, L c i uste r/ -^totai given in columns 4 and 5 of Table 7, are in the range of 0.25-0.5; the 
uncertainties are dominated by the photometry and aperture corrections of the clusters, as 
they are superposed on bright local stellar background. The light ratios are slightly larger in 
the F547M band than in the F675W band, because the clusters are 0.1-0.2 mag bluer than 
their diffuse stellar background. This color difference can be caused by an age difference 
of a few Myr, assuming that the clusters and the underlying star clouds have the same 
initial mass function. As the precise ages are unknown, we cannot model the star clouds to 
determine their masses; therefore, the cluster-to-total light ratio can be considered only as 
an approximation of the fraction of massive stars formed in clusters. In the four regions we 
analyzed, about 25-50% of the massive stars are formed in R136-class clusters. 

The Ha images in Figure 14 show that the H II regions around the clusters have rather 
well-defined boundaries where the surface brightness drops off sharply. Such morphology 
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suggests that the H II regions are likely ionization-bounded, or optically thick to ionizing 
radiation. We have measured the Ha fluxes 3 of these four regions using the apertures marked 
on the Ha images in Figure 14 and described in column 6 of Table 7. The continuum- 
subtracted Ha images are used for the flux measurements, and extinction corrections are 
made. Assuming a 10 4 K optically thick H II region, the derived Ha luminosity, Lh», can 
be used to determine the required ionizing luminosity, Q(H°), through the relation 

Q(H°) = 7.4 x 10 11 L Ua photons s _1 , (1) 

where Lh« is in units of ergs s -1 . The resultant ionizing luminosities (Qhii) of the four 
regions are given in column 7 of Table 7. 

The ionizing luminosity expected from the clusters at different ages can be calculated 
using the Starburst99 models. The ages of the clusters in the four selected regions are < 5 
Myr, as indicated by the associated bright H II regions. Unfortunately, during the first 5 
Myr a cluster's ionizing luminosity decreases rapidly, dropping from the maximum at ~ 1 
Myr to a factor of 5-7 lower at 5 Myr (Leitherer et al. 1999); therefore, the uncertainty in 
cluster age directly propagates into the uncertainty in ionizing luminosity of a cluster. We 
have adopted a cluster age of 3 Myr and calculated the expected ionizing luminosities of the 
clusters (Qciustcr) an d their ratios to those required by the surrounding H II regions. These 
results are given in columns 8 and 9 of Table 7. These ratios, 0.2-0.6, can be viewed as a 
very crude approximation of the fractions of massive stars formed in clusters. An interesting 
corollary of this result is that the ionizing luminosity, or Ha luminosity, of an H II region is 
not a sufficient diagnostic for the existence of SSCs because the majority of stars may reside 
outside clusters. 

In the four regions of star formation we considered, the fraction of massive stars in 
clusters estimated from the cluster-to-total light ratio is, within the uncertainties, consistent 
with that estimated from the cluster-to-H II-region ionizing luminosity ratio - no more than 
about half of the massive stars are formed in the R136-class clusters. On the other hand, 
some R136-class clusters, such as #4 and #5 in NGC5471D, are not superposed on a bright 
stellar background, and constitute the dominant components in their associated episode of 
star formation. The fraction of stars formed in clusters must cover a range, which varies 
according to the physical conditions of star formation. 



3 Owing to the ~ 300 km s _1 redshift of M101, two corrections need to be considered. First, the filter 
transmission of the red-shifted Ha line is ~ 93% of the peak transmission, thus the extracted Ha flux should 
be multiplied by a correction factor of 1.07. Second, the [N IIJA6548 line is red-shifted into the Ha bandpass 
at ~ 91% of the peak transmission and needs to be removed. The [N II] contamination, estimated from the 
[N II] /Ha ratios reported by Kennicutt & Garnett (1996), amounts to 1-3% of the Ha flux in most cases. 
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5.3. Interstellar Environments of the GHRs 

We examine the distribution of interstellar molecular clouds and H I gas in the three 
M101 GHRs in order to gain insight on the cluster formation process. Molecular CO ob- 
servations of these GHRs have been made by Giannakopoulou-Creighton, Fich, & Wilson 
(1999) using both single-dish telescopes and interferometers. NGC5461 has the strongest 
CO emission among the three GHRs, and its CO peaks appear concentrated toward the 
peaks of the Ha emission (see Figure 4 in Giannakopoulou-Creighton et al. 1999), where the 
massive clusters are located. The CO emission toward NGC5462 is detected in single-dish 
observations but not in interferometric observations, indicating that the molecular gas is 
distributed over a scale larger than the synthesized beam, ~3". No CO emission is detected 
in NGC5471, which may be attributed to its low metallicity. In the two GHRs where CO is 
detected, the distribution of molecular clouds is similar to that of clusters: concentrated in 
NGC5461 and distributed in NGC5462. 

The distribution of H I gas in M 101 and its relation with ionized gas have been reported 
by Smith et al. (2000, see their Figures 1-2). Their H I map shows that NGC5461 and 
NGC5462 are in the same spiral arm. Assuming a trailing arm, the offsets of H I ridge 
downstream from the stars in NGC 5461 and NGC 5462 are consistent with the expectations 
of star formation triggered by density waves (Roberts 1969). NGC 5471 has a concentration 
of H I gas but the large-scale distribution of H I show a complex inter-arm structure, which 
may have resulted from tidal interactions during the last 10 9 yr (Waller et al. 1997, and 
references therein). As H I gas can be produced by photodissociation of the natal molecular 
clouds (Allen, Atherton, & Tilanus 1985, 1986; Smith et al. 2000), converted to H II by pho- 
toionization, and dispersed by fast stellar winds and supernova explosions, the distribution 
of H I does not provide adequately pertinent information about the cluster formation. 

To study the physical conditions for cluster formation, it is necessary to examine the 
interstellar environment of the youngest clusters before the interstellar conditions have been 
altered by stellar energy feedback. The embedded young clusters that are observable in the 
infrared but not yet in the optical wavelengths (e.g., Kobulnicky & Johnson 1999; Turner, 
Beck, & Ho 2000; Johnson & Kobulnicky 2003) provide promising locations to study the 
physical conditions of cluster formation. 

5.4. Cluster Luminosity Function 

The luminosity functions (LFs) of young compact clusters have been studied in various 
types of galaxies with different star formation rates as a means to gain insight into the cluster 
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formation process. To first order, the measured LFs for young compact cluster systems in 
merging or starburst galaxies are remarkably universal, and can be approximated by a power- 
law of the form dN(L)/dL oc L a , with the exponent a ~ — 2 ± 0.2 (see Whitmore 2003, and 
references therein). The cluster LFs for a sample of nearby non-starburst spiral galaxies 
also show similar a, —2.0 to —2.4 (Larsen 2002). It has been suggested that this roughly 
universal LF is the result of fractal structure in turbulent gas (Elmegreen & Efremov 1997). 
Since the three Ml 01 GHRs show different age and spatial distributions of clusters, it is 
interesting to intercompare their cluster LFs and see if they also follow the universal cluster 
LFs. 

The LFs of clusters in NGC5461, NGC5462, and NGC5471 are presented in Figure 15. 
These LFs are constructed using raw M F 547m (without extinction correction). No complete- 
ness correction to the LFs is needed because the cutoff of our sample at the faint end, 
^F547M < 21.3, is much brighter than the detection limit, mF547M ~ 25.5. Note however that 
the two faintest bins, Mp547M = —8.0 to —9.0, should be viewed with caution, as some of 
the "clusters" may be spurious as discussed in §5.1. The number of clusters in each GHR 
is modest, so we have also constructed a combined cluster LF of the three GHRs, shown 
in the bottom panel of Figure 15. We have carried out linear least-squares fits to the log- 
arithmic LFs of clusters for the three GHRs individually and combined. The logarithmic 
LFs of clusters in the individual GHRs appear to have different slopes; the best-fit slopes for 
NGC5461, NGC5462, and NGC5471 are -1.5±0.3, -3.0±0.2, and -1.9±0.4, respectively. 
The logarithmic LF of clusters in all three GHRs has a best-fit slope of —2.3 ± 0.1. 

Compared with the universal cluster LFs, NGC5461 and NGC5471 are on the flatter 
side, and NGC5462 is on the steeper side. The small difference between NGC5461 and 
NGC5471 is not statistically significant, as each has only a small number of clusters and 
the numbers of clusters in the brightest bins are only 1-2. On the other hand, it may be 
statistically significant that NGC 5462 has a much steeper LF than NGC 5461 and NGC 5471, 
or NGC 5462 has a larger proportion of low-mass clusters. It is possible that the slope of 
a cluster LF varies according to the interstellar environment at the time when clusters 
were formed. As discussed in §5.3, the current molecular environments of NGC 5461 and 
NGC 5462 are quite different, with molecular CO highly concentrated in NGC 5461 and 
diffuse in NGC 5462. If the current environments reflect the conditions when the clusters 
were formed, the clusters in NGC 5462 would have been formed in a lower-pressure, lower- 
concentration interstellar environment. The association between a steep cluster LF and a 
low-pressure, low-concentration star formation condition is consistent with the previously 
suggested hypothesis that massive clusters are formed in high-pressure, high-concentration 
molecular clouds. 
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The cluster LF for NGC5461, NGC5462, and NGC5471 combined has an a within the 
range of the universal value, —2 to —2.4. If clusters in spiral or starburst galaxies are formed 
under conditions similar to NGC5461 and NGC5462 in random proportions, the cluster LFs 
should show a larger range of a. The scarcity of cluster LFs with a steeper than —2.5 might 
be an observational effect to some extent. Surveys of clusters in galaxies preferentially detect 
the most luminous clusters that are likely formed in high-pressure environments and follow 
an LF similar to those of NGC5461 and NGC5471; therefore, the cluster LF of a galaxy 
would be biased toward a = —2. 



5.5. Evolutionary Aspects of the Clusters 

The cluster mass, age, and size distribution of a cluster system may be used to investigate 
the dynamic evolution of clusters. Recent studies of rich clusters in the LMC have shown 
that the spread in core radius increases with cluster age, suggesting that all clusters were 
formed with small core radii but subsequently some experienced core expansion while others 
did not (Elson, Freeman, & Lauer 1989; Mackey & Gilmore 2003). It would be interesting 
to examine the clusters in M101 GHRs to see whether they follow the same core radius-age 
relation. 

The mass, age, and core radius of the LMC clusters have been derived by Mackey & 
Gilmore (2003), using the Fioc & Rocca-Volmerange (1997) population synthesis code, the 
Kroupa, Tout, & Gilmore (1993) IMF slope, and a stellar mass range of 0.1-120 M . To 
compare the M101 clusters with the LMC clusters, we have followed the Mackey & Gilmore 
(2003) method and re-estimated the masses of M101 clusters, using the same code, IMF 
slope, and stellar mass range. The new cluster mass estimates are ~ 2-3 times as high as 
the cluster masses estimated earlier in this paper, owing to the addition of stars in the mass 
range of 0.1-1 M . The M101 clusters are not sufficiently resolved for measurements of their 
core radii (r c ); therefore, we have adopted the relation r c ~ 0.35i? c g- derived from a King 
profile with a concentration parameter of c = 1.5, a median value for LMC clusters (Elson 
et al. 1987). The core radii of M101 clusters thus estimated are 0.25-1 pc, and may be 
uncertain by up to a factor of 2, if its concentration parameter spans the same range as that 
in the LMC clusters. 

To compare the M101 clusters with the LMC clusters, we present a 3-D diagram of 
cluster mass, age, and core radius in Figure 16. The data of the LMC clusters are adopted 
from Mackey & Gilmore (2003). The M101 clusters are plotted in open rhombuses, while 
the LMC clusters are plotted in filled ellipses, with the R136 cluster in a larger ellipse for 
easy identification. The M101 cluster masses were estimated using 7-pc-radius apertures, so 
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we repeated Mackey & Gilmore's derivation of the mass of R136 using this larger aperture 
and obtained a mass that is 40% higher, shown as a large open ellipse in Figure 16. While 
the LMC clusters span a large age range and show an increasing spread in core radii with 
the cluster age, the M101 clusters we analyzed are all young and small, sharing a similar 
parameter space with the R136 cluster. Among the small number of M101 clusters with 
size measurements, the younger clusters are generally smaller, but the accuracy of the size 
measurements is too limited by the linear resolution for definitive conclusions. 

Finally, we discuss the disruption time of our M101 clusters and investigate whether 
dynamical evolution has caused a significant mass change in these clusters. The disruption 
time can be assessed empirically from the break in slopes of the logarithmic age distribution 
of clusters (Boutloukos & Lamers 2003), or from comparisons with iV-body simulations 
(Baumgardt & Makino 2003). Recent studies using the empirical method have shown that 
the disruption time for a 10 4 M cluster varies greatly among galaxies, from 10 7 to 10 10 
yr, with the shortest being ~30-40 Myr in M82 and at 1-3 kpc from the nucleus of M51 
(Boutloukos & Lamers 2003; de Grijs et al. 2003). Our M101 clusters in GHRs are apparently 
younger than these disruption timescales. Thus we determine the disruption time, Tdi S , 
following the simulations of Baumgardt & Makino (2003): 



where R G is the galactocentric radius, Vq is the circular velocity in a galaxy, N is the number 
of stars in a cluster, (3 ~ 1-2, 7 = 0.02, and x ~ 0.8. For a cluster with a Salpeter mass 
function and a stellar mass range of 1-100 M , the average mass of a star is 3.09 M and 
the number of stars is iV = cluster mass/(3.09 M ). NGC5461 and NGC5462 have R G ~ 10 
kpc and V G = 185 km s -1 , and NGC 5471 has R G ~ 25 kpc and V G = 195 km s _1 (Roberts & 
Rots 1973). A 10 4 M cluster in NGC 546 1/NGC 5462 or NGC 5471 would have disruption 
times of 2.4-4.8 xlO 9 and 5.8-11.6 xlO 9 yr, respectively. The real disruption time must be 
shorter because the above estimates do not take into account processes that are important in 
disrupting clusters in GHRs, i.e., interactions with other clusters and with giant molecular 
clouds. In cases where massive clusters are concentrated in a small volume, such as the 
core of NGC 5461, cluster merger is a more important dynamic process and operates in a 
much shorter timescale than tidal disruption (Bonnell et al. 2003). Future simulations of 
dynamical evolution of clusters in GHRs using realistic conditions are needed. 
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6. Summary 

GHRs contains high concentrations of massive stars; thus, they provide an excellent 
laboratory to study modes of massive star formation and possible sites to form globular 
clusters. We have selected three very luminous but morphologically different GHRs in M101, 
NGC5461, NGC5462, and NGC5471, to determine their cluster content in order to under- 
stand cluster formation in different environments. We have obtained HST WFPC2 images 
of these GHRs with the F547M and F675W continuum filters and the F656N Ha filter. 
The continuum images are used to identify cluster candidates in each GHR and to carry 
out photometric measurements, and the Ha images are used to examine the distribution of 
interstellar gas and to determine the ionizing flux requirement. 

We have used the Starburst99 (Leitherer et al. 1999) and BC03 (Bruzual & Chariot 
2003) population synthesis models to compute the colors and magnitudes of clusters of 
different ages and masses. The colors of a cluster are dependent on its age; however, our 
two continuum passbands are not blue enough to be sensitive to young massive stars for 
age determination. Therefore, we use the distribution of ionized interstellar gas to estimate 
the approximate cluster ages, then compare the measured colors and magnitudes of cluster 
candidates to the synthetic evolutionary tracks to determine their masses. To avoid confusion 
by luminous single supergiants, only cluster candidates more luminous than M F5 4 7M = —9.0 
are analyzed for their masses. 

NGC5461 is dominated by a very luminous core, and has been suggested as a likely 
host of SSCs (Kennicutt & Chu 1988; Luridiana & Peimbert 2001). Our observations show 
that the core of NGC5461 contains three R136-class clusters superposed on a bright stellar 
background in a small region ~32 pc across. It is possible that the three R136-class clusters 
will dynamically interact and merge into an SSC. If NGC5461 were at a distance > 20 
Mpc, the clusters at its core would appear as a single cluster, and the total light would 
be Mp547M — —13.0, corresponding to a mass of ~ 10 5 M , reaching those of SSCs. It is 
possible that some of the previously reported SSCs at large distances are actually made up 
by tight groups of R136-class clusters similar to those in NGC5461. 

NGC5462 consists of numerous loosely-distributed H II regions that are individually 
much fainter than 30 Dor. Its clusters also show a loose distribution across the GHR. 
NGC5462 has the largest number of clusters among the three GHRs studied, but most of 
the clusters are older than 10 Myr and fainter than Mp547M = —9.0. 

NGC5471 contains multiple bright H II regions, some of which are comparable to 30 
Dor. A large number of cluster candidates are identified in NGC5471; the majority of the 
clusters are fainter than M F5 47 M = —9.0 and they are in bright H II regions. The mass 
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determination for clusters in NGC5471 is problematic because the observed cluster colors 
are bluer than those spanned by the synthetic cluster evolutionary tracks for Z = 0.2Z & , 
possibly as a result of uncertainties in stellar evolution models at low metallicities. The 
cluster masses are thus estimated from the magnitudes alone and may be subject to large 
errors. 

The most massive clusters in the three GHRs are in the mass range of ~ 1-3 xlO 4 
M Q , similar to R136. Two clusters in NGC5471 might be more massive as they are not 
surrounded by H II regions and are each as luminous as R136; these two may be the most 
massive clusters in the three GHRs studied. No SSCs are present in any of the three GHRs. 
We have also estimated the sizes of some clusters on their PC images, using the routine 
developed by Larsen (1999). The effective radii of these clusters are in the range of 0.7-2.9 
pc, ~ 1-3 times that of R136 (Mackey & Gilmore 2003). 

To understand the makeup of the faint cluster candidates, we have simulated WFPC2 
images of 30 Dor and NGC604 at the distance of M101. We find that single supergiants 
similar to R131 and OB associations as those in NGC604 may contribute to the population 
of faint clusters (M F5 4 7M > —9.0) in distant galaxies, while clusters similar to Hodge 301 or 
OB associations similar to LH 99 will be too faint or too extended to be identified as clusters 
even in M101. 

The three M101 GHRs show different cluster LFs. The cluster LFs of spiral galaxies can 
be described by a power-law with the exponent a in the range of —2.0 to —2.4 (Larsen 2002). 
We find that the cluster LFs of NGC 5461 and NGC 5471 are on the flatter side of the range, 
but the number of clusters is small in each GHR. NGC 5462 has the largest number of clusters 
and its cluster LF is significantly steeper, with a = —3.0 ±0.2. It is possible that the clusters 
in NGC 5462 were formed in a low-pressure, low-concentration interstellar environment. The 
combined cluster LF of the three GHRs has an a of —2.3 ± 0.1, well within the range for 
those of spiral galaxies. The universality of cluster LFs may be a statistical result from a 
cluster population with an observational bias toward the most luminous clusters. 

The distribution of molecular clouds is concentrated in NGC 5461 and diffuse in NGC 5462, 
similar to the spatial distribution of their clusters. The diffuse interstellar environment and 
the larger proportion of low-mass clusters (steep cluster LF) of NGC 5462 qualitatively sup- 
port the hypothesis that massive clusters are formed in high-pressure, high-concentration 
interstellar medium. 

We have estimated the fraction of massive stars formed in clusters using (1) clusters' 
contribution to the total stellar continuum, and (2) comparison between the ionizing flux 
expected from the clusters and the ionizing flux required by the associated H II region. Both 
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methods show that < 50% of massive stars are formed in R136-class clusters. Consequently, 
the Ha luminosity of an H II region does not provide a sufficient diagnostic for the existence 
of SSCs. 
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Fig. 1. — The POSS-II red image of M101 from the Digitized Sky Survey. The three luminous 
GHRs studied in this paper, NGC5461, NGC5462, and NGC5471, are marked. 

Fig. 2.— Color composite of HSTWFPC2 images of NGC 5461, with F547M in blue, F675W 
in green, and Ha in red. North is up and east to the left. The field-of-view is 45" x 45", or 
1.6 kpc x 1.6 kpc. 

Fig. 3.— Color composite of HSTWFPC2 images of NGC 5462, with F547M in blue, F675W 
in green, and Ha in red. North is up and east to the left. The field-of-view is 56" x 56", or 
2.0 kpc x 2.0 kpc. 

Fig. A.— Color composite of //5TWFPC2 images of NGC 5471, with F547M in blue, F675W 
in green, and Ha in red. North is up and east to the left. The field-of-view is 51" x 51", or 
1.8 kpc x 1.8 kpc. 

Fig. 5.— HST WFPC2 images of the main body of NGC 5461 in (a) F547M, (b) F675W, (c) 
Ha, (d) Ha-subtracted F675W, (e) continuum-subtracted Ha, and (f) F547M bands. Note 
that this field-of-view is smaller than that shown in Figure 2. The Ha images (c) and (e) are 
presented in different stretches to show bright and faint features, and similarly the F547M 
images (a) and (f) to show bright and faint stars/clusters. The two brightest H II regions 
from Hodge et al. (1990) are marked in (e), and cluster candidates are marked in (f). 

Fig. 6. — HST WFPC2 images of the main body of NGC 5462, displayed in a format identical 
to that of Figure 5. Note that this field-of-view is smaller than that shown in Figure 3. The 
three brightest H II regions from Hodge et al. (1990) are marked in (e), and cluster candidates 
are marked in (f). 

Fig. 7. — HST WFPC2 images of NGC 5471, displayed in a format identical to that of 
Figure 5. Note that this field-of-view is smaller than that shown in Figure 4. The five 
brightest components from Skillman (1985) are marked in (e), and cluster candidates are 
marked in (f). 
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Fig. 8. — M F 547M versus (M F 547m — -^F675w) diagram of cluster candidates in NGC5461. 
Observations of the clusters are plotted in filled circles. Additional circles are draw around 
clusters that are embedded in bright H II regions. Evolutionary tracks generated from 
Starburst99 and BC03 for a Salpeter initial mass function and a metallicity of Z = 1 Z & 
are plotted in solid and dotted curves, respectively. Ages in Myr are marked along the 
evolutionary tracks. To avoid crowding, the Starburst99 evolutionary track is shown for a 
cluster mass of 5 x 10 4 M and the BC03 track for 1 x 10 4 M . The reddening vectors of 
associated H II regions are plotted in dashed arrows, and the possible dereddening vectors 
are marked with solid arrows for clusters brighter than M F 547m = —9.0. The R136 cluster is 
plotted as a reference point: reddened R136 in a filled triangle and dereddened R136 in an 
open triangle. 



Fig. 9. — Mp547M versus (MF547M — Mf675w) diagram of cluster candidates in NGC5462. 
Symbols are the same as in Figure 8 with the addition of a dashed circle for the cluster 
candidate in an interstellar shell. The evolutionary tracks are the same as in Figure 8. 



Fig. 10. — Mp547M versus (Mp547M — A^F675W') diagram of cluster candidates in NGC5471. 
Symbols are the same as in Figure 8. The evolutionary tracks are generated from population 
synthesis models using the same parameters as those generated for NGC 5461 and NGC 5462, 
but with a metallicity of Z = 0.2 Z Q . 



Fig. 11.— (a) HST WFPC2 PC image of the clusters in the core of NGC 5461 (i.e., H1105) 
in the F547M band. The cluster numbers given in Figure 5 are again marked, (b) Binned 
HST WFPC2 WFC image in the F547M band for the same region to simulate a WFC image 
at 20 Mpc. At such a large distance, the three clusters at the core of NGC 5461 are no longer 
distinguishable from a single super-star cluster. 



Fig. 12.— HST WFPC2 PC images of the brightest clusters in NGC 5462 and NGC 5471 for 
which cluster sizes were measured: (a) NGC 5462-6, (b) NGC 5471-4 & -5, (c) NGC 5471-9, 
and (d) NGC 5471-16. The pixel size of a PC image is 070455, corresponding to a linear size 
of 1.6 pc in M101. The field-of-view of each image is 2'.'4 x 274. 
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Fig. 13.— (a) Binned MCELS green band image of 30 Dor and (b) binned HST WFPC2 
F547M images of NGC 604 to simulate WFPC2 images of 30 Dor and NGC 604 at the distance 
of M101. (c)-(e) HST WFPC2 F547M images of M101 GHRs NGC 5461, NGC 5462, and 
NGC 5471. The cluster candidates with MF547M > —9.0 are marked. All images have the 
same linear field-of-view, 350 pc x 350 pc. The four F547M images are displayed with the 
same intensity scale, while the green band image of 30 Dor is displayed with an intensity 
scale matching that of the F547M images as much as possible. 



Fig. 14. — HST WFPC2 images of four regions with R136-class clusters coexistent with 
bright diffuse stellar background in the F547M band (upper panels) and Ha band (bottom 
panels): (a) NGC 5461-1, (b) NGC 5461-8, 9, 10, (c) NGC 5462-18, and (d) NGC 5471-16. 
Circles in the F547M images mark the apertures used to measure the total stellar continuum 
emission of each region, and ellipses in the Ha images mark the apertures used to measure 
the Ha fluxes of the associated H II region. The field-of-view of each image is 8" x 8". 



Fig. 15. — M F5 4 7M luminosity functions (LFs) of candidate clusters in NGC 5461, NGC 5462, 
NGC 5471, and of the combined sample of all three GHRs. Solid lines are cluster LFs and 
dotted lines are power-law fits to these LFs. The best-fit a value and the number of cluster 
candidates are labeled in the upper left corner of each panel. 



Fig. 16. — A 3-D diagram of age, core radius, and mass of clusters in the LMC and M101. 
The M101 clusters are shown as open rhombuses, while the LMC clusters (Mackey & Gilmore 
2003) are shown as filled ellipses. The R136 cluster is shown as a larger filled ellipse and open 
ellipse for the mass estimates from Mackey & Gilmore (2003) and this study, respectively. 
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Table 1. Properties of the Three Luminous GHRs in M101 and 30 Dor in the LMC 





NGC5461 


NGC5462 


NGC 5471 


on t ~\ 

30 Dor 


Angular Size 


40" x 25" 


48" x 33" 


17" x 17" 


20' x 20' 


Linear Size a (pc) 


1400 x 875 


1680 x 1150 


600 x 600 


290 x 290 


L Ha b (ergs s" 1 ) 


2.7 x 10 40 


1.3 x 10 40 


2.2 x 10 40 


3.9 x 10 39 


Location 


in spiral arm 


in spiral arm 


outlier 


above one end of 










the LMC bar 


Ha Morphology 


one dominant core 


weak cores with 


multiple cores 


one core with 




with filaments & 


long filaments & 


with filaments 


bright loops & 




small cores around 


loops extending out 


around 


filaments around 



a We adopted the distances of 7.2 Mpc to M101 (Stetson et al. 1998) and 50 kpc to the LMC 
(Feast 1999). 

b The L Ra of NGC 5461, NGC 5462, and NGC 5471 are measured using HST WFPC2 Ha images 
in this study, and the L Hq of 30 Dor is adopted from Kennicutt & Hodge (1986). Note that these 
Lh q are not corrected for extinction. 

Table 2. Table of Observations 



Obs. Date 



Object 


(y/m/d) 


Filter 


Camera 


Exp. Time 




NGC 5461 


1999/03/24 


F547M 


WF2 


600s x2, 100s x2, 


20s xl 






F547M 


PCI 


20s x2 








F675W 


WF2 


400s x2, 50s x2, 


10s xl 




1999/03/23 


F656N 


WF2 


600s x2, 160s xl 




NGC 5462 


2000/02/01 


F547M 


WF2 


600s x2, 100s x2, 


20s xl 






F547M 


PCI 


20s x2 








F675W 


WF2 


400s x2, 50s x2, 


10s xl 






F656N 


WF2 


600s x2, 160s xl 




NGC 5471 


1997/11/01 


F547M 


WF3 


600s x2, 100s x2, 


20s xl 






F547M 


PCI 


20s x2 








F675W 


WF3 


400s x2, 50s x2, 


10s xl 






F656N 


WF3 


600s x2, 180s xl 
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Table 3. Photometry of Candidate Clusters in NGC5461 



ID «J2000 $.12000 MF547M ^F547M ~ ^F675W ^F547M ~ ^F675W 



1 


14 


03 


39 


.84 


54 


18 


56. 


.2 


-9 


.35±0. 


.01 


0. 


45±0. 


.05 


0. 


23±0 


.02 


2 


14 


03 


39 


.91 


54 


18 


56. 


.2 


-8. 


.19±0. 


.05 


1. 


26±0. 


.06 


0. 


55±0. 


.06 


3 


14 


03 


40 


.52 


54 


18 


58. 


.8 


-8. 


,97±0. 


.01 


0. 


17±0. 


.03 


0, 


16±0. 


.02 


4 


14 


03 


40 


.54 


54 


18 


59. 


.2 


-8. 


.11±0. 


.03 


0. 


14±0. 


.08 


0, 


11±0. 


.04 


5 


14 


03 


40 


.98 


54 


19 


02 


.1 


-8. 


.25±0. 


.02 


0, 


11±0. 


.06 


0, 


09±0 


.03 


6 


14 


03 


41 


.15 


54 


19 


04. 


,5 


-9 


.61±0. 


,03 


0. 


,21±0. 


.09 


0. 


10±0. 


.04 


7 


14 


03 


41 


.22 


54 


18 


57. 


.0 


-8. 


,01±0. 


.02 


0. 


17±0. 


.02 


0. 


13±0. 


.02 


ga 


14 


03 


41 


.36 


54 


19 


03 


.7 


-10 


.40±0 


,10 


0. 


59±0. 


.13 


0. 


36±0 


.10 




















-10 


.19±0. 


,06 














9 a 


14 


03 


41 


.40 


54 


19 


03 


.8 


-10 


.17±0. 


.12 


0. 


39±0. 


.16 


0. 


24±0. 


.14 




















-9 


.75±0. 


.11 














10 a 


14 


03 


41 


.42 


54 


19 


04. 


.0 


-10 


.81±0. 


,06 


0. 


27±0. 


.09 


0. 


,22±0. 


.07 




















-10 


.57±0 


.06 














11 


14 


03 


41 


.58 


54 


19 


04. 


.0 


-9 


.27±0. 


.05 


1. 


37±0. 


.05 


0. 


58±0. 


,05 


12 


14 


03 


41 


.58 


54 


19 


07. 


.8 


-8. 


,37±0. 


,01 


0. 


38±0. 


.03 


0. 


38±0 


.02 



a The photometry is given in two rows, with the first row measured with the WFC images 
and the second row measured with the PC images. 



Table 4. Photometry of Candidate Clusters in NGC 5462 



ID «J2000 £j2000 -^F547M -^F547M ~ ^F675W ^F547M _ Mf675W 



1 


14 


03 


51. 


.80 


54 


21 


52. 


.6 


-8 


.83±0. 


.04 


0. 


.32±0 


,05 





,16±0. 


,04 


2 


14 


03 


51. 


.83 


54 


21 


46. 


.2 


-8 


.08±0. 


,01 





,31±0. 


.02 





,27±0. 


.02 


3 


14 


03 


51. 


.98 


54 


21 


46. 


.7 


-8 


.64±0. 


.01 





,36±0. 


.02 





.33±0 


.02 


4 


14 


03 


52. 


.07 


54 


21 


49. 


,1 


-8 


.32±0 


,01 


0. 


,46±0. 


.01 





.43±0. 


.01 


5 


14 


03 


52. 


.38 


54 


21 


49 


,5 


-8 


.33±0 


,02 





.20±0. 


.03 





15±0. 


.03 


6 


14 


03 


52. 


.41 


54 


21 


49. 


.1 


-9 


.56±0. 


.01 





.43±0. 


.01 





.40±0. 


.01 


7 


14 


03 


52 


.42 


54 


21 


49. 


,4 


-8 


.16±0. 


,05 





,47±0. 


.05 





,44±0. 


.05 


8 


14 


03 


52. 


.43 


54 


21 


50. 


.0 


-8 


.86±0. 


.02 





.30±0. 


.02 





,28±0. 


.02 


9 


14 


03 


52. 


.84 


54 


21 


54. 


.5 


-8 


.05±0. 


,02 





,14±0. 


.03 





,12±0. 


.03 


10 


14 


03 


52. 


.86 


54 


21 


59. 


,3 


-8 


.13±0. 


.02 





.63±0 


,02 


0. 


,59±0 


,02 


11 


14 


03 


52. 


.95 


54 


21 


54. 


.2 


-8 


.55±0. 


,01 





,65±0. 


.02 





,52±0. 


.02 


12 


14 


03 


52. 


.96 


54 


22 


06 


,4 


-8 


.75±0. 


.03 


1 


,41±0. 


,03 





,37±0. 


,03 


13 


14 


03 


53 


.01 


54 


22 


00 


.6 


-8 


.44±0. 


.02 





,51±0. 


.02 





,45±0. 


.02 


14 


14 


03 


53 


.03 


54 


21 


56. 


.1 


-8 


.10±0 


,02 





,23±0 


.03 





.20±0. 


.03 


15 


14 


03 


53 


,35 


54 


22 


00 


.5 


-8 


.87±0. 


.01 





,20±0. 


.01 





17±0. 


.01 


16 


14 


03 


53 


.41 


54 


21 


59 


.1 


-8 


.61±0. 


.01 





,41±0. 


.02 





,39±0. 


.02 


17 


14 


03 


53 


.58 


54 


22 


04. 


.0 


-8 


.03±0 


.02 





.44±0. 


.03 





,42±0. 


.02 


18 


14 


03 


53 


.78 


54 


22 


11. 


.1 


-9 


.03±0 


.02 





.28±0. 


.08 





,03±0 


.03 


19 


14 


03 


53 


.98 


54 


21 


56. 


.1 


-8 


.24±0. 


.01 





,50±0. 


,02 





,48±0. 


,02 


20 


14 


03 


54. 


.00 


54 


22 


07 


,9 


-8 


.35±0. 


.02 


0. 


,19±0. 


.03 





,18±0. 


.03 


21 


14 


03 


54 


,10 


54 


22 


02 


,5 


-8 


.30±0 


.01 





,54±0. 


,02 





,41±0 


,02 


22 


14 


03 


54. 


.18 


54 


22 


06 


,6 


-8 


.24±0. 


,02 





.33±0 


.03 





,31±0. 


.03 


23 


14 


03 


54 


,19 


54 


22 


11. 


.2 


-9 


.17±0. 


.01 


0. 


,55±0. 


,01 





,51±0. 


,01 


24 


14 


03 


54 


.32 


54 


22 


09 


.1 


-8 


.09±0. 


.05 





10±0. 


.06 





,07±0. 


.06 


25 


14 


03 


54. 


.74 


54 


21 


53 


.7 


-8 


.20±0 


.01 





,27±0. 


.01 





,24±0. 


.01 
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Table 5. Photometry of Candidate Clusters in NGC 5471 



ID «J2000 <5j2000 -^F547M ^F547M — ^F675W ^F547M ~ ^F675W 



l a 


14 


04 


28.64 


54 


23 


51.9 


-8.81±0.09 


1.24±0.12 


0.41±0.11 
















-8.23±0.12 






2 a 


14 


04 


28.64 


54 


23 


52.1 


-9.31±0.08 


1.46±0.10 


0.76±0.10 
















-9.05±0.09 






3 


14 


04 


28.86 


54 


23 


48.2 


-9.08±0.03 


-0.18±0.06 


-0.13±0.04 


4 


14 


04 


28.88 


54 


23 


47.8 


-9.86±0.02 


0.15±0.03 


0.12±0.03 


o 


1 A 


04 


98 QO 


^A 


9^ 


AR A 


in 97+n ni 


_n n9+D D9 


_n fi7+n 09 


6 


14 


04 


29.10 


54 


23 


41.7 


-8.59±0.04 


0.34±0.09 


-0.09±0.05 


7 


14 


04 


29.15 


54 


23 


45.9 


-9.00±0.06 


0.04±0.09 


-0.11±0.07 


8 


14 


04 


29.17 


54 


23 


45.4 


-8.87±0.06 


1.12±0.08 


0.14±0.07 


9 


14 


04 


29.27 


54 


23 


51.2 


-9.38±0.01 


0.26±0.03 


0.19±0.01 


10 


14 


04 


29.29 


54 


23 


52.4 


-8.71±0.08 


0.99±0.11 


0.14±0.13 


11 


14 


04 


29.29 


54 


23 


52.9 


-8.21±0.05 


0.86±0.09 


0.22±0.08 


12 


14 


04 


29.33 


54 


23 


47.2 


-9.06±0.02 


0.37±0.04 


0.19±0.03 


13 


14 


04 


29.37 


54 


23 


46.5 


-8.47±0.04 


0.34±0.04 


-0.04±0.04 


14 


14 


04 


29.38 


54 


23 


46.2 


-8.08±0.06 


0.72±0.08 


0.35±0.09 


15 


14 


04 


29.39 


54 


23 


51.8 


-8.02±0.02 


0.54±0.14 


0.02±0.09 


16 


14 


04 


29.47 


54 


23 


46.4 


-10.06±0.06 


0.16±0.13 


-0.13±0.06 


17 


14 


04 


29.53 


54 


23 


46.1 


-8.68±0.14 


1.47±0.15 


0.25±0.15 


18 


14 


04 


29.54 


54 


23 


45.8 


-8.40±0.11 


1.38±0.13 


0.52±0.12 


19 


14 


04 


29.56 


54 


23 


47.5 


-8.97±0.02 


0.44±0.07 


0.21±0.03 



a The photometry is given in two rows, with the first row measured with the WFC images 
and the second row measured with the PC images. 



-36- 



Table 6. Physical Properties of Massive Clusters in M101 GHRs 





/\ge 


Su£ir niircH~QQ Ix/T^md 

kJ u til YJ LL1 o u <_/ <_/ iviaoo 


RC03 Mass 


R 




PI net or TFl 


{MjT) 


( v i n 4 m \ 

^ All) 1K/0J 


( v i n 4 m \ 


VP C ) 


iveiiiciiKS 


NGC 5461-1 


< 5 


2±0.5 


2±0.5 






NGC 5461-6 


< 5 


1±0.5 


< 1 


0.8 ±0.2 




NGC 5461-8 


< 5 


1.5±0.5 


1.5±0.5 


0.7 ±0.2 




NGC 5461-9 


< 5 


1.5±0.5 


1±0.5 






NGC 5461-10 


< 5 


3±1 


2±1 


2.1 ±0.1 


asymmetric morphology 


NGC 5461-11 


< 5 


~ 0.5 


< 0.5 






NGC 5462-6 


> 10 


~ 2 


1.5-2 


2.3 ±0.3 


asymmetric morphology 


NGC 5462-18 


< 5 


< 1 


0.5-1 






NGC 5462-23 


> 10 


~ 2 


1.5-2 






NGC 5471-2 


< 5 


~ 0.2 


~ 0.2 






NGC 5471-3 


> 10 


~ 0.5 


~ 0.5 






NGC 5471-4 


> 10 


> 1 


> 1 


2.9 ±0.3 


asymmetric morphology 


NGC 5471-5 


> 10 


~ 2 


~ 2 


1.4 ±0.1 


asymmetric morphology 


NGC 5471-7 


< 5 


- 0.5 


- 0.5 






NGC 5471-9 


> 10 


~ 1 


~ 1 


0.2 ±0.1 


probably a star 


NGC 5471-12 


< 5 


~ 0.5 


- 0.5 






NGC 5471-16 


< 5 


1.5±0.5 


1.5±0.5 


1.1 ±0.1 





Table 7. Fractional Contribution of Clusters to Total Stellar Light & Ionizing Luminosity 



Region 


Cluster 


Continuum 
Aperture 

Diameter 


-^cluster 
^total 

in F547M 


-^cluster 
^total 

in F675W' 


Ha 
Aperture 

Diameter 


Qhii 

(xlO 50 s- 1 ) 


^cluster 

(xlO 50 s- 1 ) 


a 

Qclustcr 

Qhii 


1 


NGC 5461-1 


1'.'3 


0.47 ±0.05 


0.39 ±0.04 


2'.'1 X 2'.'0 


14 


8 


0.57 


2 


NGC 5461-8,9,10 


2'/ 2 


0.43 ± 0.08 


0.38 ±0.08 


3^'5 x 2/8 


97 


24 


0.25 


3 


NGC 5462-18 


1'.'4 


0.41 ± 0.04 


0.35 ±0.04 


2'! A x 2'.' 3 


12 


4 


0.33 


4 


NGC 5471-16 


l'/5 


0.30 ± 0.03 


0.24 ±0.03 


2'.' 5 x 2'.' 3 


37 


8 


0.22 



a The ionizing luminosities of these clusters are estimated assuming a cluster age of 3 Myr. Since these clusters are only know to 
have ages < 5 Myr, the ionizing luminosity of the clusters and thus the ratio < ^g^ stcr could range from 1/4 to 1.5 times the value in 
the table. 
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